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Shell model calculations for Ni, Cu and Zn isotopes by modifying fpg interaction due 
to Sorhn et. al, [Phys. Rev. Lett. 88, 092501 (2002)] have been reported. In the present 
work 28 two body matrix elements of the earlier interaction have been modified. Present 
interaction is able to explain new experimental results for this region. 
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1. Introduction 

The nickel isotopes {Z = 28) cover three doubly-closed shells with number N — 28, 
N = 40, = 50 and therefore a unique testing ground to investigate the evolution 
of shell structure. The ^^Ni and its neighboring attracted the interest of recent re- 
search to answer the ma gicity ver sus superfluidity question related to doubly magic 
character of this nuclei. lH^EIlin] Pqj. ^t^q copper isotopes, the important question 
is related to rapid reduction in the energ y of 5/2^ state as the filling of neutrons 
started in the 1^59/2 orbital I 6 | 7 | 8 | 9 | 10 [n] ^^le Cu isotopes as 40, the state 
reveal three types of structures - single-particle, collective or coupling of single pro- 
ton with neighboring even Ni isotopes. In the Zn isotopes with two protons more 
than the semi-magic Z = 28 nickel nuclei, the spherical = 40 gap may not be 
strong enough to stabilize the nuclei in spherical shape when protons are added to 
the ^^Ni core. Thus B{E2) value is larger in the chain of Zn isotopes. To understand 
the evolution of nuclear structure, th e irnpq rtance of monopole term from the tensor 
force is pointed out by Otsuka et al. 

The experimental E{2^) for Z = 28 to Z = 36 are shown in Fig. 1. In this figure 
E{2^) states in Zn are overall lower compared to Ni and an additional decrease of 
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Fig. 1. Systematic of the experimentally observed E(2^) for Z = 28 to Z = 36 near the 
N = 28, 40 and 50 shell closure. 

the £{2^^) energy in Zn isotopes is obtained between A^=40-50 compared to 7V=28- 
40. This figure also reveal the enhancement of collectivity beyond Z = 30, because 
of rapid decrease of i?(2+). Below the Ni chain, the spectroscopy of Cr, Fe isotopes 
shown an incr eased collectivity toward = 40, revealing the collapse of this shell 
closure. H^Ell 

Following our recent shel l-mode l (SM) studies for neutron-rich F isotopes 
odd and even isotopes of Fe odd-odd Mn isotopes odd-mass 6i,63,65(2;q 

isotopes^^ni^ and odd-even ''^"^^Ga isotopesl^, in the present work, large scale shell 
model calculations have been performed for neutron rich Ni, Cu and Zn isotopes for 
40 < < 50 in fpggf2 model space. The low-lying energy levels and B{E2) values 
have been calculated and compared with the recent experimental data. 

The paper is organized as follows: In Section 2 details of calculation are de- 
scribed. In Section 3 results and discussion are presented. Finally in Section 4 we 
give conclusions. 

2. Details of Calculation 

Large scale shell model calculations have been performed for neutron rich nickel, 
copper and zinc isotopes with 40<N<50 using '^''Ca as a core by including the /7/2 
orbit. In the previous work^^Sl^ large scale shell model calculations were performed 
for neutron rich nickel, cop per and zinc isotopes for 40< N <50, using three differ- 
ent versions | 23 | 24 | 25 | 26 | 27 | effective interactions for the model space consisting of 
the P3/2, Pi/2, /5/2 and 59/2 . Both, however, yield unsatisfactory results in certain 
aspects, viz. 

(i) large E(2+) value for very neutron rich nuclei (^^Ni and ^"Zn), 

(ii) small B{E2) values in comparison to experimental values, and 
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Shell model calculations in Ni region 3 

(iii) for ^^Cu, ""^Cu and ^^Cu, the ground state is 3/2^ as compared to the experi- 
mental indication of 5/2^. 

This is probably due to neglect of the fi/2 orbit for the protons. In view of 
this, /7/2 orbit has been included in the valence space of protons by taking ^"Ca 
as core. Effective interaction for fpg valence space for both protons and neutrons 
with ^°Ca as core has been constructed by Sorlin et al.^ The main drawback of 
this interaction is that the effective proton single-particle energies of f'j/2 orbital 
become lower than /5/2 which is not realistic. In the present work, we have modified 
relevant matrix elements to account for this discrepancy. 



2.1. Model space 

In the present work, we have used fpg model space comprising of the O/7/2, lp3/2j 
0/5/2, lpi/2 active proton orbitals and O/7/2, 1^3/2, O/5/2, lpi/2, 0^9/2 neutron 
orbitals with eight /7/2 frozen neutrons, more accurately '^^Ca core with eight f7/2 
frozen neutrons. The single-particle energies for O/7/2, IP3/2, O/5/2, lpi/2, O59/2 
orbitals are 0.0, 2.0, 6.5, 4.0 and 9.0 MeV respectively. The relative single-particle 
energies are taken from the excitation energies of the low-lying negative parity states 
in 49Ca.l2] 



2.2. Effective Interaction 

For the fpg valence space, an effective interaction with *°Ca as core has been 
reported in Ref.^ This interaction has been built using fp two-body matrix elements 
(TBME) from Ref. ISH! and 1^3/2, O/5/2, lpi/2, O59/2 TBME from Ref. ^4 For the 
common active orbitals in these subspaces, matrix elements were taken from. ^211 
As the latter interaction has been defined for a ^^Ni core, a scaling factor of A~^/^ 
amplitude is applied to take into account the change of radius between the ^°Ca and 
^^Ni cores. As more and more neutrons are added in the /7/2 shell, the excitation 
energy of the 9/2+ state decreases as reflected in Fig. 2 for ''^^''^Sc isotopes. The 
remaining matrix elements are taken from 77/259/2 TBME from Ref. 



2.3. Monopole correction to the single-particle levels 

As more and more neutrons are added in gg/2 orbital, /7/2 orbital becomes more 
repulsive and /5/2 more attractive. Thus proton f^/2 is pulled down while /7/2 is 
lifted up, as N increases, this is due to monopole interaction produced by tensor 
force between a proton in j>.< = I ±1/2 and a neutron in j> < = ^ '±1/2. Effects 
on P3/2 and pi/2 orbitals are small and can be neglected. In the present work, we 
have modified 59/2/7/2 matrix elements by subtracting 100 keV and 59/2/5/2 matrix 
elements by adding 100 keV. We have also check other set of values like 50 and 150 
keV to modify 59/2/7/2 and 59/2/5/2 set of matrix elements. The results with 100 
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Fig. 2. Effective single-particle energies in Sc isotopes using pfg9 interaction. 

keV modification is more reasonable. This modified interaction is named as pfg9a. 
The effective single-particle energies for ^^~^^Sc isotopes are shown in Fig. 3. 



2.4. State-of-the-art in the calculations 

The modified interaction has been modified to reproduce the experimental data of 
Cu isotopes. These being odd nuclei, the available experimental data is sparse. The 
aim of this tuning is to first reproduce the ground state properties of ''^"''^Cu and 
then apply it to Ni and Zn isotopes. In the original fpg interaction, level crosses 
the fj/2 ^°Sc and become lower in energy. This appears to be unrealistic. 

In this modified interaction the /7/2 effective single-particle energy is always 
lower than which is expected. Since the dimension of matrices involved for 
pfg9/2 space is very large, truncation of the full shell- model space are necessary. 
We have allowed neutron: t^, jumps from (P3/2, f5/2i P1/2) to 59/2, protons: i^r 
jumps from /7/2 to (^3/2, /5/2, P1/2) orbital for each nucleus and each isotopes, 
an example of this truncation are shown in Fig. 4. We performed calculations for 
maximum dimension which are feasible by SGI-Cluster computer at GANIL. Up to 
this dimension, the states seems to be more convergent, which can be taken as the 
final result. 

All the SM calculations have been carried out using the code ANTOINE I ^4 | 30 | 31 | 
at SGI-Cluster computer at GANIL and DGCTIC-UNAM computational facility 
KanBalam. 
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Fig. 3. Effective single-particle energies in ^^Sc isotopes using pfg9a interaction. In this figure, 
the results for pfg9 interaction are also shown as dotted lines. 
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Fig. 4. Model space and truncation, neutron: t^ jumps from (p3/2:/5/2iPl/2) to gg/2, protons: 
U jumps from to (P3/2>/5/2,Pl/2) orbital. 

2.5. Binding Energies 

To compare our shell model results with the experimental binding energies relative 
to binding energies of '^'^Ca we use following formula: 



February 22, 2012 1:11 WSPC/INSTRUCTION FILE ws-mpla 



6 P.C. Srivastava 



Table 1. The experimental (Expt.) and shell model (SM) g.s. energies. The values of energies are 
in MeV. 



'^«Ni40 


™Ni42 


72Ni44 


'*Ni46 


76Ni48 




Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 




-248.3 -227.7 


-260.2 -239.5 


-271.1 -249.9 


-281.7 -259.5 


-293.7 -267.3 






"CU42 


''■^CU44 


''='CU46 


''''CU48 




Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


-257.9 -235.1 


-271.0 -248.2 


-283.4 -259.7 


-294.7 -270.4 


-305.4 -279.4 


-315.2 -287.2 


'^"Zn40 


^2Zn42 


'^4Zn44 


■'''Zn46 


''«Zn48 


™Zn50 


Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


Expt. SM 


-269.0 -242.6 


-283.7 -257.4 


-297.5 -271.8 


-310.0 -283.6 


-321.3 -293.5 


-332.0 -302.4 







Eb = -B, = E{SM) +Ec + Em. (1) 
In this Coulomb energies relative to "'"Ca given by^^ 

7r(7r — 1) .1 , 

Ec^e^TT + V^TT ^ h Ki/TTiy + [-ttJ&c, (2) 

where 7r(i/) stands for the number of valence protons (neutrons) for ^°Ca core. We 
have taken the values of parameters Ctt, V^tttt, K-i^ and be from Refl^^l These param- 
eters are determined by fitting to the measured Coulomb displacement energies of 
Ni(Z=28)< Z < Mo(Z=42) with 32< N < 50. The resulting parameters are 

= 9.504 MeV, = 0.228 MeV, V^^ = -0.036 MeV, be = 0.030 MeV. (3) 

The monopole expression is given by 1331 

Em = e^n + a^n{n - 1) + b{T{T + 1) - ^n). (4) 

where n is total number of valence particles and T is the total isospin. The is an 
average particle core interaction and a and b are the isoscalar and isovector global 
monopole correction. The values of the parameters at A = 42 are ^321 

= -861 ± 0.01 MeV, a = 0.041 ± 0.003 MeV, b = 0.119 ± 0.006 MeV. (5) 

In Table 1, we have tabulated experimental (Expt.) and shell model (SM) g.s. 
energies. 



3. Results and Discussion 

The yrast levels for the nickel, copper and zinc isotopes for 40< N < 50 are shown 
in Figs. 5, 6 and 7. 
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Fig. 5. Experimental 1351 and calculated ground-state energies for Cu isotopes from A= 



to 79. 



3.1. Copper isotopes 

For ^^Cu, the ground state 3/2" is well predicted by pfgQa interaction. The cal- 
culated 1/2" and 5/2" state is reverse in comparison to experimental levels. The 
calculated 1/2" and 5/2" state is 134 keV and 378 keV lower in energy from the 
experimental levels. The order of 7/2^, 7/2^ and 9/2" is well reproduced by this 
interaction. The 7/22^ is only 62 keV below than experimental value. In ^^Cu, the 
3/2" is interpreted by 7rp3/2, the 5/2" is interpreted by 7rp5/2- The first 7/2" is 
interpreted by 7rp3/2®2"''(^*Ni) and the second 7/2" is interpreted by t^p^^2- 
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Fig. 6. Yrast levels of ^^Cu isotopes with pfg9a interaction. 
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The experimental and calculated ground-state energies for Cu isotopes are shown 
in Fig. 5. In Fig. 6, the yrast levels of ^^~^^Cu isotopes with pfg9a interaction are 
shown. For ^^Cu, the ground state 3/2~ is well predicted by pfg9a interaction. The 
order of 5/2^ and l/2~ state is reverse in comparison to experimental levels. The 
calculated 1/2^ level is 303 keV higher than the experimental value, while 5/2~ 
level is 179 keV lower than the experimental result. The order of 7/2^ , 7 j%2. ^'^^ 
9/2~ is well reproduced and are higher in energy from the experimental values. In 
the '^"'^Cu, the 3/2 ~ is interpreted by 7rp3/2, and the 5/2" is interpreted by 
The first 7/2" is interpreted by 7rp3/202+(''°Ni) and the second 7/2" is interpreted 
by Trp-j^. 

For ^'^Cu, pfg9a interaction predicts 5/2" state lower in energy by 15 keV from 
3/2" state, while experimentally 3/2" is ground state. 1/2" state is 462 keV higher 
in energy from experimental value. The 9/2" lies in between two 7/2" states. In 
the ^'^Cu, the 3/2" is interpreted by Trp3/2, and the 5/2" is interpreted by 'Kp^/2- 
The first 7/2" is interpreted by 7rp3/2(8)2+('^^Ni) and the second 7/2" is interpreted 
by Trp-j^. 

Recently, there is experimental indication of 5/2" as a ground state in ^^Cu at 
REX-ISOLDE, CERN. iSlThe pfg9a interaction also predicts 5/2" as a ground 
state. For '"^Cu and ™Cu, this interaction predicts 5/2" as a ground state which is 
also expected from experiment. The calculations predict that the 1/2", 9/2" and 
7/2" states are at very high in energy for ''^Cu isotopes. 

3.2. Nickel isotopes 

The nickel isotopes (Z = 28) cover three doubly-closed shells and therefore, a unique 
testing ground for large scale shell model calculations. Experimentally, the 8+ iso- 
merism though expected to be present in a whole chain from ^^Ni to ''^Ni was found 
to be suddenly absent in ^^Ni and ''^Ni and present in ''°Ni and ''^Ni. The ground 
state spin from ^^Ni to ^®Ni is correctly reproduced by pfg9a interaction. 

The experimental and calculated ground-state energies for Ni isotopes are shown 
in Fig. 7. In Fig. 8, the yrast levels of ^^"'^^Ni isotopes with pfg9a interaction are 
shown. For ™Ni, the first excited state 2+ is calculated at 1.516 MeV which is 
257 keV above than the experimental level. The states 4+ at 2.229 MeV, 6+ at 
2.678 MeV, and 8+ at 2.860 MeV are well predicted at 2.677 MeV, 2.958 MeV 
and 3.095 MeV. The agreement with the experimental data is reasonable. For ^^Ni 
the first excited state 2+ is calculated at 1.362 MeV which is 266 keV above the 
experimental level. The states 4+ at 1.941 MeV and 6+ at 2.396 MeV, are well 
predicted at 2.248 MeV and 2.470 MeV. Theoretically, there is 8+ state at 2.621 
MeV, but experimentally their is no indication of 8+ isomeric states. For ^*Ni, only 
two excited states are experimentally known. The first excited state 2+ is calculated 
at 1.534 MeV which is 510 keV above the experimental value. The state 4"'" at 1.763 
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Fig. 7. Experimental and calculated ground-state energies for Ni isotopes from A=68 to 76. 
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Fig. 8. Yrast levels of °° '°Ni isotopes with pfg9a interaction. 



MeV is predicted at 2.354 MeV. The calculated 6+, 8+ state is at 2.612 MeV and 
2.740 MeV respectively. For ^^Ni, the calculated first 2+ state is slightly higher in 
energy about 557 keV from the experimental value. The states 4+ at 1.922 MeV, 
6+ at 2.276 MeV and isomeric states 8+ at 2.420 MeV are well predicted at 2.169 
MeV, 2.420 MeV and 2.532 MeV. The agreement for 4+, 6+ and 8+ states with the 
experiment is good. 

3.3. Zinc isotopes 

The experimental and calculated ground-state energies for Ni isotopes are shown 
in Fig. 9. In Fig. 10, the yrast levels of ^°~*°Zn isotopes with pfg9a interaction 
are shown. The ground state spin for ^'^"^"Zn isotopes for 40<N<50 is correctly 
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Fig. 9. Experimental 1331 and calculated ground-state energies for Zn isotopes from A=70 to 80. 

reproduced by pfg9a interaction. 

For ^"Zn, the first excited 2+ state is calculated at 0.603 MeV which is 282 keV below 
than experimental value. The states 4+ at 1.79 MeV, 6+ at 2.89 MeV and 8+ at 3.75 
MeV are well predicted at 1.47 MeV, 2.37 MeV and 2.53 MeV. The calculated values 
are compressed in comparison to the experimental value. For ^^Zn, experimentally 
only 2+ state at 0.653 MeV is known, which is predicted well by pfg9a interaction 
with energy difference of 74 keV. The calculated values of states 4+, 6+ and 8+ are 
at 1.178, 1.923 and 2.601 MeV respectively. For '^"Zn, the first excited 2+ state is 
calculated at 0.825 MeV which is 219 keV above the experimental value. The state 
4+ at 1.419 MeV is predicted at 1.551 MeV. The states 6+ and 8+ are at 2.783 
MeV and 3.430 MeV. For ™Zn, the first excited 2+ state is calculated at 0.997 MeV 
which is 398 keV above the experimental value. The first 4+ state is calculated with 
energy difference of 679 keV from the experimental value. The calculated 4+, 6+ 
and 8+ states are at 1.975 MeV, 3.036 MeV and 3.243 MeV respectively. For '^^Zn, 
the first excited 2+ state is calculated at 1.307 MeV, which is 477 keV above the 
experimental value. The states 4+ at 1.621 MeV, 6+ at 2.528 MeV and 8+ at 2.673 
MeV are predicted at 2.247 MeV, 2.923 MeV and 2.955 MeV. Experimentally, the 
difference between 6"*" and 8"'" states is 145 keV while theoretically, it is only 32 keV. 
For ''®Zn, the first excited 2+ state is calculated at 1.953 MeV which is 461 keV 
smaller than the experimental value. 

The calculated and experimental E(2j'") and E(4^) for Ni and Zn isotopes are shown 
in Fig. 11. The high values of E(2^) for Ni at N=40 is an indication of shell closure at 
N=40. For Zn isotopes, the E(2f ) and E(4^) at N=40 and 50 are high in comparison 
to neighboring isotopes reflects shell closure at N=40 and N=50. 

3.4. The B{E2) systematics in the copper isotopes 

The calculated values of first excited states of 1/2^, 5/2^ and 7/2~ with its cor- 
responding experimental values are shown in Fig. 12 for ^^^^^Cu isotopes. In the 
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Fig. 10. Yrast levels of ™ '"Zn isotopes with pfg9a interaction. 
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Fig. 11. The calculated and experimental E{2^) and E(4^) for Ni and Zn isotopes as a function 
of neutron number. 



lower part of this figure, the calculated and experimental B{E2) values correspond- 
ing to different transitions namely, B{E2; 1/2" 3/2"), B{E2; 5/2" 3/2") and 
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B{E2; 7/2^ 3/2^) are also shown. The calculated and experimental B{E2) val- 
ues are tabulated in Table 2. The previously calculated value by Stefanescu et al. 
are shown within bracket.^ The calculated B{E2; 1/2^ — > 3/2^) values show much 
better agreement with experimental data compared to those of Stefanescu. The low 
B{E2) value beyond N=40 for 5/2~ state confirms its 7r/5/2 single-particle charac- 
ter. The sharp drop in the excitation energy of 5/2~ state beyond N=40 could be 
due to the monopole migration. The large i3(i52)value for l/2~ state depart it from 
the single-particle character of 7rpi/2 type. The small B{E2;7 /2~ — > 3/2~) values 
for 73,75,77,79^^ isotopes can be seen from wavefunctions of 7/2^ and 3/2^ states in 
Table 3, which show single-particle character of these states. In the first part of the 
Fig. 12, for the E{l/2~) state the difference between predicted and experimental 
energy for ^^Cu, ''-^Cu, '''^Cu isotopes increases and its corresponding B{E2) value 
also increases. 





Fig. 12. The calculated and experimental B{E2) in W.u. for ^^Cu isotopes. The shell model 
B(E2) values were calculated with the standard effective charges e,r = 1.5e and = 0.5e. 



3.5. The B{E2) systematics in the nickel isotopes 

The B{E2;2^ — > 0^) value in ^^Ni is much lower as compared to ^^Ni. This low 
value was interpreted by Sorlin et al. ^ as originating from the enhanced neutron 
pair scattering at iV = 40, which is referred to as superfluid behavior of the neutrons. 
The constancy of the B{E2) values beyond iV = 40 can be understood as follows: 
E2 strength as np — nh excitation across = 40 for odd n cannot contribute due to 
parity conservation and for even n are dominated by pair scattering. Langanke et al. 
'^performed microscopic calculations of the B(E2) in even-even nickel isotopes and 
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Fig. 13. The calculated and experimental B(E2;2'[ — > Oj ) in W.u. for Ni and Zn isotopes as 
a function of neutron number. The shell model B(E2) values were calculated with the standard 
effective charges e-,^ = 1.5e and = 0.5e. 



found that the small observed B{E2) value is not necessarily an argument for a 
shell closure at N=40, but it simply reflects the fact that the lowest 2+ state in ^^Ni 
is primarily a neutron excitation. Van de Walle et al.^^showed that B{E2) values 
in the Ni chain show a parabolic evolution between two magic numbers (N=28 and 
40), hinting a seniority-like behavior. Perm et al. ^ measured high B{E2) value 
for ^°Ni at GANIL, and attributed it due to rapid proton core polarization when 
neutrons are added to the lg9/2 orbit. The calculated B{E2;2f 0^) values for 
68-76]\jj j^j-g shown in Fig. 13, the B{E2) value is almost constant from 68-74]\jj .^^^^ 
its value again decreases for ^^Ni. This decrease in B{E2) value is probably due to 
oncoming of the next shell closure at N=50 for ''^Ni. 

3.6. The B{E2) systematics in the zinc isotopes 

The experimental B{E2;2^ 0^) values in the Zn isotopic chain show a similar 
trend towards iV = 40 as the Ni isotopes up to ^®Zn. Though for ™Zn ( at = 
40), the B{E2) value suddenly increases. Leenhardt et al.^ give three effects for 
supporting this increased collectivity: the addition of two protons out of the Ni core, 
the maximum in neutron pairing correlations at = 40, and the presence of the 
strongly downslopping I = 4 Nilsson neutron orbitals close to the Fermi surface. 
Kenn et a/. ^3 indicate that the inclusion of the l(?g/2 orbit in the valence space is 
important in order to reproduce the increased B{E2; 2^ 0^) values in the ™Zn. 

The calculated B{E2; 2^ 0^) values from present work using pfg9a interac- 
tion for ^°~^°Zn are shown in Fig. 13. The calculated value show similar trends as 
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experimental values from ™Zn to ^°Zn. The B{E2) value at N=50 is very low, this 
is an indication of shell closure at this value. 

Table 2. B{E2) values (in W.u.) for ^^"^^Cu isotopes. The calculated results shown with corresponding value in the bracket from earlier work of Stefanescu 
The shell model B{E2) values were calculated with the standard effective charges = 1.5e and ei, = 0.5e. 

''^CU40 ^^CU42 ^3CU44 ■^5CU46 ^^Cu48 '^^Cuso 

Expt. pfg9a Expt. pfg9a Expt. pfg9a Expt. pfg9a Expt. pfg9a Expt. pfg9a 

B{E2;l/2- ^3/2-) 10.4[1.0] 11.5(7.1) 20.4[2.2] 12.5(7.3) 23.1[2.1] 12.41(7.5) 13.9 7.5 9.2 

B(£;2;5/2- ^. 3/2") 3.0[0.3] 9.9(1.6) 3.9[0.5] 2.7(1.7) 4.4[0.5] 2.07(1.3) 2.7 2.1 1.4 

B(£2;7/2- ^ 3/2-) 4.6[0.7] 7.9(1.2) 10.7[1.2] 1.2(1.5) 14.9[1.8] 0.0038(2.3) 0.0032 0.0002 0.0016 



Table 3. The calculated wave- function components of the lowest 5/2 , 3/2 , 7/2 , and 1/2 states in Cu isotopes. 
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^(/7/2/5/2)^(/7/2P3/2/5/2Pl/259/2) 


62% 
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(^/7/2/5/2)^{/7/2P3/2/5/2Pl/259/2) 
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■''■(/7/2Pl/2)^(/7/2P3/2/5/2Pl/2ff9/2) 
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4. Conclusion 

In the present work, the existing pfg% interaction for fpgg/2 space with ^"^Ca core 
has been modified by changing 71/5/2*^39/2 and 71/7/21^39/2 matrix elements. The new 
interaction is named as pfg9a has been tuned for Cu isotopes and tested for Ni and 
Zn isotopes. The pfg9a interaction give 0+ ground state in an even-even nucleus 
like in Ni and Zn isotopes which is the characteristic of any reasonable interaction. 
These results indicate that further modification in interaction and even inclusion 
of 1^5/2 orbit is important in the shell model calculations. The modification of 
Sorlin et. ai, interaction is also attempted in,ESIi3ut it will be remain to test this 
interaction is universal or not for this region. An attempt by including 1^5/2 orbit 
in fpgg /2 space to explain collectivity for fp shell nuclei is recently reported in. 1221 
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